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Recent Progress on Tau Lepton Physics 
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Some important aspects of hadronic r decays are reviewed: the determination of a 3 from the inclusive r hadronic 
width, the measurement of \V US \ through the Cabibbo-suppressed decays of the r, and the theoretical description 
of the t — > Vt-Ktt spectrum. The present status of other relevant electroweak topics, such as charged-current 
universality tests or bounds on lepton-flavour violation, has been already summarized in ref. [1]. 



1. The inclusive hadronic width of the tau 

The hadronic r decays turn out to be a beau- 
tiful laboratory for studying strong interaction 
effects at low energies [2, 3]. The r is the 
only known lepton massive enough to decay into 
hadrons. Its semileptonic decays are then ideally 
suited for studying the hadronic weak currents. 

The inclusive character of the total r hadronic 
width renders possible an accurate calculation of 
the ratio [4-8] 



R T 



r[r- 



v T hadrons] 



r[r- 



= Rt Y + Rt a + R t 



The theoretical analysis involves the two-point 
correlation functions for the vector Y$ = 'ipjj^ipi 
and axial-vector = ipj^^fsipi colour-singlet 
quark currents (i, j = u, d, s): 



d*xe^(0\T(jMx)J?A0n0),(l 



which have the Lorentz decompositions 



IT 



ij,J 



(?) = (-^V + ^)n« (q 2 ) 



(2) 



where the superscript (J — 0, 1) denotes the an- 
gular momentum in the hadronic rest frame. 

The imaginary parts of IIi/j-(<7 2 ) are propor- 
tional to the spectral functions for hadrons with 
the corresponding quantum numbers. The semi- 
hadronic decay rate of the r can be written as 
an integral of these spectral functions over the 



invariant mass s of the final-state hadrons: 



R T = l2ir l _ 



ds 



in 



m- 



1 + 2 At ) Imn (1) (s) + Imn (0) (s) 



(3) 



The appropriate combinations of correlators are 
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The contributions coming from the first two terms 
correspond to R T ,v and R t ,a respectively, while 
Rt,s contains the remaining Cabibbo-suppressed 
contributions. 

The integrand in Eq. ([3]) cannot be calculated 
at present from QCD. Nevertheless the integral it- 
self can be calculated systematically by exploiting 
the analytic properties of the correlators n( J )(s). 
They are analytic functions of s except along the 
positive real s-axis, where their imaginary parts 
have discontinuities. R T can then be written as a 
contour integral in the complex s-plane running 
counter-clockwise around the circle \s\ = m 2 .: 



R T = &m (t> I 1 — - 
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This expression requires the correlators only 
for complex s of order m 2 ., which is signifi- 
cantly larger than the scale associated with non- 
perturbative effects. Using the Operator Product 



1 



2 



Expansion (OPE), U^(s) = J2d C d /(~ s ) D/2 > 
to evaluate the contour integral, R T can be ex- 
pressed as an expansion in powers of 1/m 2 .. The 
uncertainties associated with the use of the OPE 
near the time-like axis are heavily suppressed by 
the presence in ([5]) of a double zero at s — rn 2 .. 
The combination R T y+A can be written as [6] 

Rt,V+A = N C \V ud \ 2 Sew {1 + S F + S NF } , (6) 

where Nq — 3 is the number of quark colours and 
Sew = 1-0201 ± 0.0003 contains the electroweak 
radiative corrections [9-11]. The dominant cor- 
rection (~ 20%) is the perturbative QCD contri- 
bution S P , which is already known to 0{a i s ) [6,12] 
and includes a resummation of the most impor- 
tant higher-order effects [7,13]. 

Non-perturbative contributions are suppressed 
by six powers of the r mass [6] and, therefore, 
are very small. Their numerical size has been 
determined from the invariant-mass distribution 
of the final hadrons in r decay, through the study 
of weighted integrals [14], 

r l ds (i-j-\ k C_lV ^l, (7 ) 

Jo \ m l) \ m %) ds 

which can be calculated theoretically in the same 
way as R T . The predicted suppression [6] of the 
non-perturbative corrections has been confirmed 
by ALEPH [15], CLEO [16] and OPAL [17]. The 
most recent analysis [18] gives 

6 NP = -0.0059 ±0.0014. (8) 

The QCD prediction for R T y + A is then com- 
pletely dominated by dp; non-perturbative effects 
being smaller than the perturbative uncertainties 
from uncalculated higher-order corrections. The 
result turns out to be very sensitive to the value 
of a s (m 2 ), allowing for an accurate determination 
of the fundamental QCD coupling [5,6]. The ex- 
perimental measurement R T y + A = 3.479 ±0.011 
implies [18] 

a s {m 2 T ) = 0.344 ± 0.005 oxp ± 0.007 th . (9) 

The strong coupling measured at the r mass 
scale is significantly larger than the values ob- 
tained at higher energies. From the hadronic de- 
cays of the Z, one gets a s (M§) = 0.1191±0.0027 
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Figure 1. Measured values of a s at different 
scales. The curves show the energy dependence 
predicted by QCD, using a s (m 2 ) as input. The 
corresponding extrapolated a s (M^) values are 
shown at the bottom, where the shaded band dis- 
plays the r decay result within errors [18]. 



[12,18,19], which differs from a s (m 2 ) by more 
than 20(T. After evolution up to the scale Mz 
[20] , the strong coupling constant in (J9j) decreases 
to [18] 

a s {M 2 z ) = 0.1212 ±0.0011, (10) 

in excellent agreement with the direct measure- 
ments at the Z peak and with a better accuracy. 
The comparison of these two determinations of a s 
in two very different energy regimes, m T and Mz, 
provides a beautiful test of the predicted running 
of the QCD coupling; i.e., a very significant ex- 
perimental verification of asymptotic freedom. 

2. Perturbative contribution to R T 

The recent calculation of the O(af) contribu- 
tion to n( 0+1 ) (s) [12] has triggered a renewed the- 
oretical interest on R T [12,18,21,22]. The pertur- 
bative contribution Sp is extracted from the Adler 



3 



function 

n=0 v 7 

For three flavours, the known coefficients take the 
values: K = K x = 1; # 2 = 1.63982; K 3 (MS) = 
6.37101 and tf 4 (MS) = 49.07570 [12]. 

The perturbative component of R T is given by 

5 P = Y J K nA (n) (a s ), (12) 

71=1 

where the functions [7] 
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are contour integrals in the complex plane, which 
only depend on a T = a s (m 2 )/ir. Using the 
exact solution (up to unknown /3 n >4 contribu- 
tions) for a s (s) given by the renormalization- 
group /3-function equation, they can be numer- 
ically computed with a very high accuracy [7]. 
One can easily check that the results are very 
stable under changes of the renormalization scale 
and rather insensitive to the truncation of the (3 
function (putting f}± = has a negligible impact). 
Thus, the resulting theoretical uncertainty on 5p 
is small. 

However if, instead of adopting the known val- 
ues for A^(a s ), one expands a s (—s) in powers of 
a 8 (m T ) inside the the integrals (TIB")) , the large log- 
arithmic running along the circle s — m 2 exp (i(j>) 
(</>e[0, 2tt]) gives rise to a nearly divergent series 
of the form Sp = J2 n =i(Kn + 9n) a ri where the 
g n coefficients depend on K m<n and on /3 m <„: 

<5 (0) = a T + 5.20 a 2 T + 26 A a 3 + 127 a\ + ■ ■ ■ (14) 

The "running" g n contributions are much larger 
than the original K n coefficients containing the 
Adler function dynamics (<?2 = 3.563, 33 = 19.99, 
#4 = 78.00) [7]. These generates a sizeable renor- 
malization scale dependence, which is much larger 
than the naively expected 0{a\) effect. The ra- 
dius of convergence of this expansion is actually 
quite small. A numerical analysis of the series [7] 



shows that, at the three- loop level, an upper esti- 
mate for the convergence radius is a TiConv < 0.11, 
which is very close to the physical value. Thus, 
the fixed-order expansion (TT4"1) should not be used 
for accurate predictions of R T . The result ([9]) 
has been correctly obtained using Eq. (fT2|) with 
the exact values of the functions A^ n \a s ). The 
slightly different results quoted in refs. [12, 21] 
originate in their use of the pathological fixed- 
order expansion p^]) Fl 

3. \V US \ determination from tau decays 

The separate measurement of the |A5| = 
and |AiS| = 1 tau decay widths provides a very 
clean determination of V us [23,24]. To a first ap- 
proximation the Cabibbo mixing can be directly 
obtained from experimental measurements, with- 
out any theoretical input. Neglecting the small 
SU(3)-breaking corrections from the m s — 
quark-mass difference, one gets: 

|K S | SU(3) = \V ud \ (^-) 1/2 - 0.210±0.003. 
\iW,V+AJ 

We have used \V ud \ = 0.97418 ± 0.00027 [25], 
R T = 3.640 ± 0.010 and the value R T . S = 
0.1617 ± 0.0040 [24], which results from the most 
recent BaBar [26] and Belle [27] measurements 
of Cabibbo-suppressed tau decays [28]. The new 
branching ratios measured by BaBar and Belle 
are all smaller than the previous world aver- 
ages, which translates into a smaller value of 
i? Tj 5 and \V US \. For comparison, the previous 
value R T , S = 0.1686 ± 0.0047 [18] resulted in 
|Ks| su(3) = 0.215 ±0.003. 

This rather remarkable determination is only 
slightly shifted by the small SU(3)-breaking con- 

1 A better convergence of the fixed-order expansion J 141 
is enforced in Ref. [21] through an artificial cancelation 
of the K n and g n contributions at higher orders. Since 
R T does not get corrections from D = 4 terms in the 
OPE, this behaviour is trivially accomplished assuming 
that the perturbative series is dominated by an n = 2 IR, 
renormalon. While this provides an interesting academic 
model of higher-order contributions, the resulting wild be- 
haviour of the Adler series is totally ad-hoc and gener- 
ates problems for weighted distributions of the form (0. 
The non-perturbative correction in JHJ would no longer be 
valid within this model, making the low value of a B (m T ) 
claimed in [21] unjustified. 
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tributions induced by the strange quark mass. 
These corrections can be estimated through a 
QCD analysis of the differences [23, 24, 29-36] 



SR 



T>kl 



IK. 1 2 



(15) 



The only non-zero contributions are proportional 
to the mass-squared difference to 2 — m d or to 
vacuum expectation values of SU(3)-breaking op- 
erators such as SO4 = {0\m s ss — mddd\0) ss 
(-1.4 ± 0.4) • 1CT 3 GeV 4 [23,29]. The dimen- 
sions of these operators are compensated by cor- 
responding powers of m 2 , which implies a strong 
suppression of 5R^} [29]: 



SIC 



hi 



24 S, 



EW 



m 2 (m 2 ) 



(1-eS) A kl (a s 



-2tt 



2 so. 



(16) 



where e d = m d /m s = 0.053 ± 0.002 [37]. The 
perturbative corrections Aki{a s ) and Qki(cx s ) are 
known to O{o 3 s ) and 0(a 2 ), respectively [29,36]. 

The J — contribution to Aoo(a s ) shows a 
rather pathological behaviour, with clear signs of 
being a non-convergent perturbative series. For- 
tunately, the corresponding longitudinal contri- 
bution to SR T = 8R®° can be estimated phe- 
nomenologically with a much better accuracy, 
6R T \ L = 0.1544 ± 0.0037 [23,38], because it is 
dominated by far by the well-known r — ► v T ir 
and r — > i/ T A' contributions. To estimate the re- 
maining transverse component, one needs an in- 
put value for the strange quark mass. Taking the 
range m s {m T ) = (100 ± 10) MeV [m s (2 GeV) = 
(96 ± 10) MeV], which includes the most recent 
determinations of m s from QCD sum rules and 
lattice QCD [38], one gets finally 5R T . th = 0.216± 
0.016, which implies [24] 



\V US 



Rt,, 



Rt,V+A 



0.2165 ± 0.0026 



th 



cxp 



1/2 



±0.0005 



th 



(17) 



A larger central value, \V US \ = 0.2212 ± 0.0031, is 
obtained with the old world average for R Tt s- 

Sizeable changes on the experimental determi- 
nation of R T ,s are to be expected from the full 



analysis of the huge BaBar and Belle data sam- 
ples. In particular, the high-multiplicity decay 
modes are not well known at present. Thus, the 
result (|17p could easily fluctuate in the near fu- 
ture. However, it is important to realize that the 
final error of the V us determination from r de- 
cay is completely dominated by the experimen- 
tal uncertainties. If R T ,s is measured with a 1% 
precision, the resulting V us uncertainty will get 
reduced to around 0.6%, i.e. ±0.0013, making r 
decay the best source of information about V us . 

An accurate measurement of the invariant- 
mass distribution of the final hadrons could make 
possible a simultaneous determination of V us and 
the strange quark mass, through a correlated 
analysis of several weighted differences SRf . The 
extraction of m s suffers from theoretical uncer- 
tainties related to the convergence of the pertur- 
bative series Aki(a s ), which makes necessary a 
better understanding of these corrections. 

4. r — > v T Kn and K — > irlDi 

The decays r — > v T Kn probe the same hadronic 
form factors investigated in processes, but 
they are sensitive to a much broader range of 
invariant masses. A theoretical understanding 
of the form factors can be achieved, using an- 
alyticity, unitarity and some general properties 
of QCD, such as chiral symmetry and the short- 
distance asymptotic behaviour [2,3]. 

Figure [2] compares the resulting theoretical de- 
scription of the t decay spectrum [39] with the 
recent Belle measurement [27]. At low values of s 
there is clear evidence of the scalar contribution, 
which was predicted previously using a careful 
analysis of Kir scattering data [38,40]. From the 
measured r spectrum one obtains Mk* — 895.3 ± 
0.2 MeV and T K . = 47.5±0.4 MeV [39]. Since the 
absolute normalization is fixed by K13 data to be 
|K s |/+°M0) = 0.21664 ±0.00048 [41], one gets 
then a theoretical prediction for the branching 
fraction, Br(r" -> v t K s tx-) = 0.427 ± 0.024%, 
in good agreement with the Belle measurement 
0.404 ± 0.013%, although slightly larger. 

The t determination of the vector form fac- 
tor f^(s) [39,42] provides precise values for its 
slope and curvature, = (25.2 ± 0.3) • 10~ 3 and 
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X'l = (12.9 ± 0.3) ■ 1(T 4 [39], in agreement but 
more precise than the corresponding Ki 3 mea- 
surements [41]. 




Vs [GeV] 



Figure 2. Theoretical description [39] (solid line) 
of the Belle t~ -> v t K s -k~ data [27]. The K* 
(dashed-dotted) and scalar (dotted) contributions 
are also shown. 
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